Quantized Quasi-Two Dimensional Bose-Einstein Condensates with Spatially 

Modulated Nonlinearity 
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We investigate the localized nonlinear matter waves of the quasi-two dimensional Bose-Einstein 
condensates with spatially modulated nonlinearity in harmonic potential. It is shown that the whole 
Bose-Einstein condensates, similar to the linear harmonic oscillator, can have an arbitrary number of 
localized nonlinear matter waves with discrete energies, which are mathematically exact orthogonal 
solutions of the Gross-Pitaevskii equation. Their novel properties are determined by the principle 
quantum number n and secondary quantum number I: the parity of the matter wave functions 
and the corresponding energy levels depend only on n, and the numbers of density packets for each 
quantum state depend on both n and I which describe the topological properties of the atom packets. 
We also give an experimental protocol to observe these novel phenomena in future experiments. 
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Introduction. — Since the remarkable experimental re- 
alization [H, [1, 0] of Bose-Einstein condensations (BEC), 
there has been an explosion of the experimental and the- 
oretical activity devoted to the physics of dilute ultra- 
cold bosonic gases. It is known that the properties of 
BEC including their shape, collective nonlinear excita- 
tions are determined by the sign and magnitude of the s- 
wave scattering length. A prominent way to adjust scat- 
tering length is to tune an external magnetic field in the 
vicinity of a Feshbach resonance [4(. Alternatively, one 
can use a Feshbach resonance induced by optical or elec- 
tric field [Bj]. Since all quantities of interest in the BEC 
crucially depend on scattering length, a tunable interac- 
tion suggests very interesting studies of the many-body 
behavior of condensate systems. 

In the past years, techniques for adjusting the scatter- 
ing length globally have been crucial to many experimen- 
tal achievements @, Q . More recently, condensates with 
a spatially modulated nonlinearity by manipulating scat- 
tering length locally have been proposed [1,0, Jj], II, 12]. 
This is experimentally feasible due to the flexible and 
precise control of the scattering length with tunable in- 
teractions. The spatial dependence of scattering length 
can be implemented by a spatially inhomogeneous exter- 
nal magnetic field in the vicinity of a Feshbach resonance 

0. 



However, so far, the studies of BEC with spatially mod- 
ulated nonlinearity are limited in the quasi-one dimen- 



sional cases [H, 0, [To], 11, [Hj. Moreover, in the study 
of nonlinear problems no one discusses their quantum 
properties which are common in linear systems such as 
the linear harmonic oscillator. In this Letter, we extend 
the similarity transformation [12[ to the quasi-two dimen- 
sional (quasi-2D) BEC with spatially modulated nonlin- 
earity in harmonic potential, and find a family of sta- 
ble localized nonlinear matter wave solutions. Similar 



to the linear harmonic oscillator, we discover that the 
whole BEC can be quantized which is unexpected be- 
fore. Their quantum and topological properties can be 
simply described by two quantum numbers. We also for- 
mulate an experimental procedure for the realization of 
these novel phenomena in 7 Li condensate [H, Q3]. This 
opens the door to the investigation of new matter waves 
in the high dimensional BEC with spatially modulated 
nonlinearities. 

Model and exact localized solutions. — The system con- 
sidered here is a BEC confined in a harmonic trap V(r) = 
m(oj\r 2 +uSlz 2 )/2, where m is atomic mass, r 2 = x 2 + y 2 , 
and co± , co z are the confinement frequencies in the radial 
and axial directions, respectively. In the mean- field the- 
ory, the BEC system at low temperature is described by 
the Gross-Pitaevskii (GP) equation in three dimensions. 
If the trap is pancake-shaped, i.e. uo z ^> cu±, it is reason- 
able to reduce the GP equation for the condensate wave 
function to a quasi-2D equation TH, [H, [l7| 



-^ xx ^ yy )^ 2 (x 2 ^y 2 )^^g(x,yM 2 ^^ (1) 



where u = uj±/cj z , the length, time and wave function 
i/j are measured in units of ah = y/h/mujz^u;^ 1 , a^ 1 
and g(x,y) = 47ra s (x,y) represents the strength of in- 
teratomic interaction characterized by the s-wave scat- 
tering length a s (x,y), which can be spatially inhomoge- 
neous by magnetically tuning the Feshbach resonances 

Now we consider the spatially localized stationary so- 
lution ip(x,y,t) = (j)(x,y)e~' lfJjt of Eq. (pQ) with (f>(x,y) 
being a real function for hm^^i^^ y) = 0. This 
maps Eq. (pp) onto a stationary nonlinear Schrodinger 
equation \<j> X x + \<l>yy- \u 2 (x 2 +y 2 )(l)-g(x, y)^ 3 ^^ = 
[lq . Here fi is the real chemical potential. Solving this 
stationary equation by similarity transformation [12[ , we 
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FIG. 1: (color online). The interaction parameter g(x,y) for 
two secondary quantum numbers: (a) I = and (b) I = 1 with 
uo = 0.02, v = 0.1. It is seen that g(x,y) is a smooth function 
when I — and develops singularity when / gets large. 



obtain a families of exact localized nonlinear wave solu- 
tions for Eq. (pQ) as 



(n + l)K(k)ri (a s _ iilt 
tp n = ' cn(0, k)e \ 



0,2,4,-. 



(2) 



FIG. 2: (color online). The density distributions of the quasi- 
2D BEC with spatially modulated nonlinearities in harmonic 
potential, for different principle quantum numbers n in Eqs. 
(2)- (3) with Eq. (4), where the secondary quantum number 
/ = 0, the parameters uo, v are 0.02 and 0.1, respectively. The 
unit of space length x,y is 1.69 fim. Figs. 2(a)-2(c) show 
the density profiles of the even parity wave function (2) for 
n = 0, 2 and 4, respectively. Figs. 2(d)-2(f) demonstrate the 
density profiles of the odd parity wave function (3) for n — 1, 3 
and 5, respectively. 



^ = (n + l)K(k)r, ^ fc)e _^ n = 1; 3; 5; . . . (3) 

where k = a/2/2 is the modulus of elliptic function, v 
is a positive real constant, K(k) = J 2 [1 — k 2 sin 2 s]d$ is 
elliptic integral of the first kind, sd = sn/dn with sn, cn 
and dn being Jacobi elliptic functions, 77 and g are de- 
termined by 

= (n + l)K(k)erf [y/fa (x + y) /2], 
77 = e w ^KummerU[-/i/(2u;), 1/2, cj (x - ?/) 2 /2], (4) 

^(x,^) = -2^^/(V)e-" ( ^ )2 , 

here erf(x) = e~ T<2 dr is error function, and 

KummerU(a, c, 5) [19] is Kummer function of the second 
kind which is a solution of ordinary differential equation 
sA (s) + (c — s)A (s) — aA(s) = 0. It is easy to see that 
when I x |, 1 2/ 1 — > 00 we have ^ n — > for solutions tjj n in 
Eqs. (2)- (3) with Eq. (4), thus they are localized bound 
state solutions. 

In the above construction, it is observed that the num- 
ber of zero points of function 77 in Eq. (4) is equal to 
that of function KummerU[— fi/(2uo), 1/2, uo (x — y) 2 /2], 
which strongly depends on uo and the ratio ji/uo. We as- 
sume the number of zero points in 77 along line y = — x is 



/. In the following, we will see that integer n is associated 
with the energy levels of the atoms and integers n, / de- 
termine the topological properties of atom packets, so n 
and / are named the principal quantum number and sec- 
ondary quantum number in quantum mechanics. In ad- 
dition, the three free parameters lj, fi and v are positive, 
so the dimensionless interaction function g(x,y) is neg- 
ative, which indicates an attractive interaction between 
atoms. There are known atomic gases with attractive in- 
teractions realized by modulating magnetic M| technique, 
for examples, the 85 Rb and 7 Li atoms |2l fl4j. 

To translate our results into units relevant to the ex- 
periments [3, [13 ]. we take the 7 Li condensate contain- 
ing 10 3 - 10 5 atoms in a pancake-shaped trap with ra- 
dial frequency uo± = 2tt x 10 Hz and axial frequency 
uo z = 2tt x 500 Hz [2l|. In this case, the ratio of trap 
frequency uo in Eq. (1) is 0.02 which is determined 
by uo±/uo z . The unit of length is 1.69 /im, the unit of 
time is 0.32 ms and the unit of chemical potential is 
nK. The spatially inhomogeneous interaction parameter 
g(x,y) is independent of principal quantum number n but 
is strongly related to the secondary quantum number 
In the Fig. 1, we show that for uo = 0.02, v = 0.1, func- 
tion g(x,y) is smooth in space when / = and develops 
singularity when the / gets large. 

Quantized quasi-2D BEC. — In order to investigate the 
quantum and topological properties of the localized non- 
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FIG. 3: (color online). The density distributions of the quasi- 
2D BEC in harmonic potential for different secondary quan- 
tum number /. Figs. 3(a)-3(d) show the density profiles of 
the even parity wave function (2) for principle quantum num- 
ber n = 0, and Figs. 3(e)-3(h) show the density profiles of 
the odd parity wave function (3) for n = 1, corresponding to 
/ = 0, 1, 2,3. The other parameters are the same as that of 
Fig. 2. 



linear matter waves in quasi-2D BEC described by Eqs. 
(2)- (3) with Eq. (4), we plot their density distributions 
by manipulating the principal quantum number n or sec- 
ondary quantum number /. 

Firstly, when the secondary quantum number / is fixed, 
we can modulate the principal quantum number n to an- 
alyze the novel matter waves in quasi-2D BEC. Fig. 2 
shows the density profiles in quasi-2D BEC with spatially 
modulated nonlinearities in harmonic potential for / = 0. 
It is easy to see that the matter wave functions in Eq. 
(2) satisfy ijj n (—x, —y) = ^ n (x, y), so they are even parity 
and are invariant under space inversion. Figs. 2(a)-2(c) 
demonstrate the density profiles of the even parity wave 
functions (2) with Eq. (4) for n = 0,2,4, which corre- 
spond to a low energy state and two highly excited states. 
The numbers of atoms N n = j f dxdy\i/j n (x : y, t)\ 2 for the 
three states are 7V = 3.76 x 10 3 ,7V 2 = 6.84 x 10 4 ,7V 4 = 
2.633 x 10 5 , respectively. The matter wave functions in 
Eq. (3) satisfy tjj n (—x,—y) = —ip n (x,y), which denotes 
that they are odd parity. Figs. 2(d)-2(f) demonstrate the 
density profiles of the odd parity wave functions (3) with 
Eq. (4) for n = 1,3,5, which correspond to three highly 
excited states. The numbers of atoms for the three states 
are N\ = 4.016 x 10 4 , N 3 = 2.493 x 10 5 , N 5 = 7.28 x 10 5 , 
respectively. It is observed that when the secondary 
quantum number / = 0, the number of nodes along line 
y = x for each quantum state is equal to the correspond- 



ing principal quantum number n, i.e. the nth level quan- 
tum state has n nodes along y = x. And the number 
of density packets increases one by one along line y = x 
when the n increases. This is similar to the quantum 
properties in the linear harmonic oscillator. 

Secondly, when the principal quantum number n is 
fixed, we can tune the secondary quantum number / 
to observe the novel quantum phenomenon in quasi-2D 
BEC. In Fig. 3 we demonstrate the density distribu- 
tions of quasi-2D BEC in harmonic potential for differ- 
ent secondary quantum number. Figs. 3(a)-3(d) show 
the density profiles of the even parity wave function (2) 
with Eq. (4) for n = 0, and / = 0,1,2 and 3, respec- 
tively. It is seen that the number of nodes for the density 
packets along line y = — x is equal to the corresponding 
secondary quantum number / which describes the topo- 
logical patterns of the atom packets, and the number 
of density packets increases one by one when / increases. 
Figs. 3(e)-3(h) show the density profiles of the odd parity 
wave function (3) with Eq. (4) for n = 1 and / = 0, 1, 2, 3. 
We see that the number of density packets increases pair 
by pair when / increases. The number of density packets 
for each quantum state is equal to (n + 1) x (7 + 1), and 
all the density packets are symmetrical with respect to 
lines y = ±x, as shown in Figs. 2-3. 

Normalization energy vs chemical potential. Next 
we calculate the normalization energy of each quantum 
states numerically. The total energy of the quasi-2D 
BEC is E{$) = J Jdxdy[\V^\ 2 + \uo 2 (x 2 + y 2 )\^\ 2 + 
\g{x, y)\^W So the normalized energy is given by 
EW)/N = fi - ^JJdxdyg(x : yM 4 with N = 
J J dxdy\^\ 2 . Fig. 4 shows the relations of the normal- 
ization energy E(ip)/N with chemistry potential for dif- 
ferent principle quantum numbers n. It is observed that 
for the fixed n, the normalization energy is approxima- 
tively linear increase with respect to chemistry potential, 
i.e., d(E(ijj)/N)/d/jL > 0. Fig. 4(a) demonstrates that the 
normalization energy for the even parity wave function 
(2) increases when the principal quantum number n in- 
creases. So does the odd parity wave function (3), as 
shown in Fig. 4(b). It is shown that the energy levels of 
the atoms are only associated with the principle quantum 
number n. These are similar to energy level distribution 
of the energy eigenvalue problem for the linear harmonic 
oscillator described by linear Schrodinger equation. 

Stability analysis. — Stability of exact solutions with re- 
spect to perturbation is very important, because only 
stable localized nonlinear matter waves are promis- 
ing for experimental observations and physical appli- 
cations. To study the stability of our exact solutions 
(2)- (3) with Eq. (4), we consider a perturbed solution 
^!/,*) = [^(^J/j+t^y,^* of Eq. (1). Here 
<f) n (x,y) are the exact solutions of the stationary non- 
linear Schrodinger equation \(j) xx + \(j> yy — \oo 2 (x 2 + 
y 2 )(j) — g(x,y)(j) 3 + fief) = 0. ^(x,y,t) <C 1 is a small 
perturbation to the exact solutions and ^(x,y,t) = 
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FIG. 4: (color online). The normalization energy E(ip)/N vs 
chemical potential /x, with N = f f dxdy\ip\ 2 . (a) even parity 
wave function (2) with principal quantum numbers n = 0, 2, 4 
and (b) odd parity wave function (3) with n — 1,3,5. Here 
the parameters cj = 0.02 and v — 0.1. 



FIG. 5: (color online). Eigenvalue for different principal quan- 
tum number n with parameters oo = 0.02, \i = 0.001 and 
v — 0.1. It is shown that only for n — 0, 1,2,3,4,5 are the 
localized nonlinear matter wave solutions (2)- (3) with Eq. (4) 
linear stability. 



[R(x,y) + I(x 1 y))e lXt is decomposed into its real and 
imaginary parts [23]. Substituting this perturbed so- 
lution to the quasi-2D GP equation (1) and neglect- 
ing the higher-order terms in (R, /), we obtain a stan- 
dard eigenvalue problem L+R = Ai", L_I = \R, where 
A is eigenvalue, R, I are eigenfunctions with L + = 
-\{dl + d 2 ) +3g(x,y)(j) n (x,y) 2 + \u?(x 2 + y 2 ) - \i and 
L_ = -\{d 2 x + d*)+g(x, y)$ n {x, y) 2 + \uo 2 (x 2 + ?/ 2 ) - /i. 
Numerical experiments show that when uj = 0.02 and /i, ^ 
are arbitrary no n- negative constants, only for principle 
quantum number n = 0, 1,2, 3, 4, 5 are the eigenvalues A 
of this eigenvalue problem real. This suggests that for 
ou = 0.02 the exact localized nonlinear matter wave solu- 
tion (2) is linear stability only for n = 0, 2, 4 and solution 
(3) is linear stability only for n = 1,3,5, see Fig. 5. It 
is seen that when the frequencies of pancake-shaped trap 
is fixed, the stability of the exact solutions (2)- (3) with 
Eq. (4) rests only on the principle quantum number n. 

Experimental protocol — We now provide an experi- 
mental protocol for creating the quasi-2D localized non- 
linear matter waves. To do so, we take the attractive 7 Li 
condensate [H, Q3], containing about 10 3 ~ 10 5 atoms, 
confined in a pancake-shaped trap with radial frequency 
oj± = 2-7T x 10 Hz and axial frequency uj z = 2tt x 500 
Hz (2l|. This trap can be determined by combination 



of spectroscopic observations, direct magnetic field mea- 
surement, and the observe d sp atial cylindrical symmetry 
of the trapped atom cloud [2l| . The next step is to realize 
the spatial variation of the scattering length. Near the 
Feshbach resonance 0, [H, [H, Hi], the scattering length 
a s (B) varies dispersively as a function of magnetic field 
B, i.e. a s (B) = a[l + A/(B - B)}, with a being the 
asymptotic value of the scattering length far from the 



resonance, Bq being the resonant value of the magnetic 
field, and A being the width of the resonance. For the 
magnetic field in z direction with gradient a along x-y 
direction, we have B = [Bq + aBi(x,y)]e* z . In this case, 
the scattering length is dependent on x and y. In real ex- 
periments, the spatially dependent magnetic field may be 
generated by a microfabricated ferromagnetic structure 
integrated on an atom chip 25|, [26[, such that interac- 
tion in Fig. 1 can be realized. In order to observe the 
density distributions in Figs. 2-3 clearly in experiment, 
the 7 Li atoms should be evaporatively cooled to low tem- 
peratures, say in the range of 10 to 100 nK. After the 
interaction parameter in Fig. 1(a) is realized by modulat- 
ing magnetic field properly, the density distributions in 
Fig. 2 can be observed for different numbers of atoms by 
evaporative cooling, for example, the numbers of atoms 
in Fig. 2(a)-2(c) are 3.76 x 10 3 ,6.84 x 10 4 , 2.633 x 10 5 , 
respectively. The density distributions in Fig. 3 can also 
be observed by changing the scattering lengths through 
magnetic field for various atom numbers. 

Conclusion. — In summary, we have discovered a new 
family of stable exact localized nonlinear matter wave 
solutions of the quasi-2D BEC with spatially modulated 
nonlinearities in harmonic potential. Similar to the linear 
harmonic oscillator, we introduce two classes of quantum 
numbers: the principle quantum number n and secondary 
quantum number I. The matter wave functions have even 
parity for the even principle quantum number and odd 
parity for the odd one, the energy levels of the atoms are 
only associated with the principle quantum number, and 
the number of density packets for each quantum state is 
equal to (n + 1) x (/ + 1). We also provide an experimen- 
tal scheme to observe these novel phenomena in future 



5 



experiments. Our results are of particular significance to 
matter wave management in high dimensional BEC. 
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